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Abstract. Using a recursive solution of the Yang-Mills equation, we calculate analytic 
expressions for the gluon fields created in ultra-relativistic heavy ion collisions at small times 
t. We have worked out explicit solutions for the fields and the energy momentum tensor up 
to 4th order in an expansion in r. We generalize the McLerran- Venugopalan model to allow 
for a systematic treatment of averaged charge densities fj 2 that vary as a function of transverse 
coordinates. This allows us to calculate radial, elliptic and directed flow of gluon fields. Our 
results can serve as initial conditions for hydrodynamic simulations of nuclear collisions that 
include initial flow. 



1. Introduction 

Nuclear collisions at high energies at the Relativistic Heavy Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) can create and probe novel forms of QCD matter. One example is a 
deconfined phase of QCD called Quark Gluon Plasma (QGP). A thermalized QGP is created in 
nuclear collisions at RHIC and LHC after a time of about 1 fm/c. Relativistic hydrodynamics 
is the tool of choice to describe the evolution of QGP in this stage [U12]. Large uncertainties in 
the hydrodynamic evolution come from the incomplete understanding of the initial conditions, 
i.e. the evolution of the nuclear collision prior to 1 fm/c. This is where another new phase of 
QCD, Color Glass Condensates (CGC) plays an important part. The natural limit of a nucleus 
at very high energy is characterized by a saturated gluon density in which the gluon field can 
be described by a quasi-classical field [3JH]. The earliest phase of the collision of two nuclei can 
thus be described as the collision of two sheets of color glass jjj]. The classical field resulting 
from the collision then decays and equilibrates to a plasma of quarks and gluons. 

To be more precise, in the color glass formalism partons are distinguished by momentum 
fraction. Sources have large Bjorken-x while the small- x gluons are treated through the classical 
field F^ u . They are related by the Yang- Mills equation [D fJi ,F flu ] = J". Before the collision, 
the two incoming nuclei on the positive and negative light cone are dominated by transverse 
chromo-electric and chromo-magnetic fields which are perpendicular to each other [3j 0] . After 
the collision, the fields in the forward light cone can be parameterized as = ±x^A and a 
transverse field A 1 , i = 1, 2. A and A % depend on the proper time r = \Jt 2 — z 2 and transverse 
coordinates but not on the space-time rapidity r\. For small r, a series expansion A^ = ^2 n T n A^ 
in powers of r can be employed. A few lowest orders in time should give an adequate and analytic 
description of the gluon fields immediately after the collision [6]. 



The collision initially (r — > ) creates longitudinal chromo-electric and magnetic fields 

E = ig[A\,A 2 ], B = -ige^[A\,A{] 

between the nuclei, where A\ and A\ are the fields of nucleus 1 or 2 before the collision in light 
cone gauge. The next order in r gives us linearly growing transverse electric and magnetic fields 

Efa = -^smhrilD^Eo] + cosh V , B }) , Bfa = ± (cosh r, e« [D^,E Q } - sinh V [D\ B }) . 

We have calculated analytic expressions for the fields up to fourth order in r explicitly and 
will discuss those results elsewhere [7J. The behavior we find is qualitatively consistent with 
numerical studies of classical field QCD [8]. 

The energy momentum tensor of the gluon field after the collision can be easily calculated. 
Up to order 0(r 2 ) — we skip higher orders here for brevity — it can be written in the form 

Bi B 2 \ 

A + D E B'Jr 

E A-D B'Jr { ' 

B'Jr B'Jr (-A + C)/r 2 J 

in the r, x, y, r\ coordinate system. Here A = {E 2 + Bq)/2 is the initial energy density at r = 0, 
and Bi and B[ are the coefficients of the Poynting vector T° l = Bi cosh r\ + B[ sinh r\ and linear 
in r [i = 1,2 are transverse indices). C, D, E tx r 2 are coefficients quadratic in r. All these 
coefficients have been calculated analytically [Tj. One can check that this energy momentum 
tensor is boost invariant owing to the original boost-invariant setup of the colliding nuclei. 
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2. Flow in the McLerran Venugopalan Model 

In the McLerran- Venugopalan (MV) model [3j[l], the simplest implementation of color glass, 
an observable O calculated from sources in the two nuclei has to be averaged over all possible 
configurations of color source distributions p. Each collision samples arbitrary p due to the short 
time scale of the collision compared to internal time scales of the nuclear wave function. The 
distribution of charge is usually assumed to be Gaussian around zero with an average value p? 
of the variation. We use the definition 

(pa(x±)p b (y±)) = -^^jP 2 (x±)5 2 (x_ L - yj . (2) 

We have suppressed the dependence on the longitudinal coordinate for brevity. Note that we 
have generalized the original expression from the MV model such that p 2 can be a function 
of transverse coordinate x±, instead of being an average charge density that is assumed to be 
homogeneous. Of course the latter is unrealistic and the generalization needs to be made if long- 
range transverse dynamics is to be described. We have shown explicitly that if p 2 is varying 
slowly such that |^ 2 (xj_)| 3> m _1 |VV 2 (^±)l ^ m~ 2 \ V l V J /i 2 (^j_)|, color glass dynamics and 
long-distance dynamics which are not described by CGC can be safely separated. Here m is an 
infrared regulator. We find that the generalized two gluon correlation function is [TJ 

(4(x m^)) = ^£0^ - jSn? + wir^" 4 " - hr^ ■ (3) 



where the first term is the same as in the original MV model [S]. Q here is the saturation scale. 
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Figure 1. Left: The hydro-like flow component B (arrows) and the underlying energy density. 
Middle: The rapidity-odd flow component B'. Right: The Poynting vector T° l at r\ = 1 together 
with the energy density. The contour lines indicate the density profile of one of the nuclei. 



Let us now focus on the transverse Poynting vector T° l in the energy momentum tensor which 
describes the energy flow of the gluon field. The term even in r\ after averaging reads 

q 6 iV c (iV c 2 - 1) , , , , Q 2 
Bi ~ ~ T 32^ V (/il/i2) ln (TiT^F ' { ) 

This expression gives radial and elliptic flow of energy oc — V*e. This is similar to what a 
hydrodynamic flow field would look like for the same energy density e if the hydrodynamic 
pressure is a monotonous function of e. The term odd in rj is 

Interestingly this term can describe a directed flow of energy. Note that this term does not 
violate boost-invariance of the energy momentum tensor. On the other hand boost-invariant 
hydrodynamics would never be able to create such a flow field. Rather one can trace this term 
back to the QCD analogon of Gauss' Law. This makes this rapidity-odd flow field a curious 
phenomenon. The flow vector fields B and B' from two gold nuclei with Woods-Saxon profiles 
/j, 2 colliding with impact parameter b = 8 fm are shown in Fig. [TJ The net Poynting vector T° l 
exhibits a clear dipole asymmetry in the transverse plane when rj ^ 0. Flow in color glass at 
midrapidity has recently also been studied numerically in |10| . 



3. Hydrodynamic Initial Conditions 

Since the classical gluon field only describes the very earliest stage of a nuclear collision we 
need to translate the energy momentum tensor into initial conditions for hydrodynamics. The 
dynamics of equilibration is outside the purview of either hydrodynamics nor classical field 
theory and would require a more complete description of quantum non-abelian dynamics. One 
can however derive a matching hydrodynamic energy momentum tensor simply using conserved 
currents both on the classical field side and the hydrodynamic side. For ideal hydrodynamics 
this was worked out in [1 1|. [P2] by just imposing energy and momentum conservation. The flow 
of energy generally translates into a flow of particles on the hydrodynamic side. The direction 
of the hydrodynamic transverse flow velocity is parallel to the Poynting vector and its size is a 
monotonous function of the magnitude of the Poynting vector. A similar correspondence can be 
worked out in the case of matching to viscous hydrodynamics [7j. 



4. Directed Flow Phenomenology 

Without a further hydrodynamic evolution we can nevertheless check wether our somewhat 
surprising direct flow result is qualitatively consistent with experiments. The first thing to 
check is the direction of the flow. It turns out that our directed flow vector B' points away from 
the spectators for a given side ?7>Oor?7<Oofa collision system with finite impact parameter. 
This is consistent with experiments. The rapidity dependence of particle directed flow v% is also 
roughly consistent with a sinhy shape as a function of particle rapidity y [TB| [Tl"]. 

Note that B' vanishes for b = if the two nuclei are equal. Interestingly, if an asymmetric 
system like Cu+Au collides with 6 = the flow vector B' becomes radial, but points outward 
in the direction of the smaller nucleus and inward in the direction of the larger nucleus. An 
asymmetric system with b ^ would exhibit a striking blend of both effects. Looking for such 
a strong dependence of flow on rapidity in central collisions of asymmetric systems could be a 
signature for the eminence of strong gluon fields early in nuclear collisions. 



5. Conclusion 

To summarize, if the McLerran Venugopalan model is generalized to allow for variations of the 
average charge density across the nucleus flow phenomena in early stages of nuclear collisions can 
be studied. Here we use analytic results for early gluon fields and their energy momentum tensor 
that can be obtained in an expansion of the Yang-Mills equations around the time of collision. 
We find radial and elliptic flow patterns similar to hydrodynamic behavior. In addition we also 
obtain a rapidity-odd term in the flow of gluon energy which resembles directed flow. We have 
argued that this flow will translate into directed flow of particles in the further course of the 
collision. The results are in qualitative agreement with results from the STAR experiment. 
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